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Abstract 



The first and second-order supersymmetry transformations are used to generate Ha- 
miltonians with known spectra departing from the trigonometric Poschl-Teller poten- 
cy ' tials. The several possibilities of manipulating the initial spectrum aie fully explored, 
and it is shown how to modify one or two levels, or even to leave the spectrum unaf- 



fected. The behavior of the new potentials at the boundaries of the domain is studied. 

> 1 Introduction 

o 

. There is a growing interest nowadays in the design of systems whose Hamiltonians have pre- 
^ ! scribed energy spectra, and the simplest technique to achieve this goal is the supersymmetric 
^ I quantum mechanics (SUSY QM) [1]. In this procedure, jieparting from an initial solvable 
^ ■ Hamiltonian H it can be constructed a new solvable one H with slightly modified spectrum, 
O ■ by using a finite-order differential intertwining operator [2-27]. The ingredients to imple- 
ment these transformations are seed solutions of the initial stationary Schrodinger equation 
^ . associated to factorization energies which do not coincide in general with the eigenvalues 
I of H. By iterating appropriately this method as many times as needed, one could construct 
Hamiltonians whose spectra are arbitrarily close to any desired one. 

In the case that the intertwining operator is of first order the procedure can be implemented 
by using as seed one Schrodinger solution which factorization energy is less than or equal 
to the ground state energy of H [2-16]. In order to surpass successfully this restriction, one 
needs to use interwining operators at least of second order [17-27]. The resulting second- 
order SUSY QM offers several interesting possibilities of spectral manipulation [15, 18, 19]: 

(i) two new levels can be placed between a pair of neighbor physical ones -Ej-i, Ei of H; 

(ii) one new energy can be created at an arbitrary position; (iii) one level can be moved; (iv) 
there is not modification of the initial spectrum; (v) one physical energy can be deleted; (vi) 
two neighbor physical levels can be deleted. 

The SUSY techniques have been extensively applied to several interesting examples for 
which the x-domain is the full real line (e.g. the harmonic oscillator) or the positive semi- 
axis (e.g. the radial oscillator or the Coulomb problem). In order to complete the scheme, 
it is important to apply them to cases where the x-domain is a finite interval, let us say 
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[xi, Xr]. An example of this kind, to be explored in detail in this paper, is the trigonometric 
Poschl-Teller potential [16,27-29]. This is closely related to several potentials widely used 
in molecular and solid state physics [28]. Since the SUSY transformations modify slightly 
the initial spectrum, it turns out that a lot of new potentials are available to be used as model 
in physical applications. 

In the next section we will survey quickly the A;-th order SUSY QM, with special emphasis 
placed in the first and second-order cases [15]. In section 3 we will build up the first and 
second-order SUSY partners of the trigonometric Poschl-Teller potential. In section 4 we 
will finish the paper with our conclusions. 



2 Supersymmetric quantum mechanics 

The study of systems ruled by the supersymmetry algebra with two generators, 

[Qi,H,,] = 0, {Qi,Qj} = 6ijH,„ t,3 = l,2, (1) 

realized in the way 

o - ^ + ^' o - ^'-^ m 

^=(50)' = (0 ^ ) ' = (^0^ 55t) ' 

where is a k-ih. order differential operator intertwining two Schrodinger Hamiltonians 

Has 



HB^ = B^H, (4) 

1 ~ 1 rf2 



H = ~— + Vix), H=—— + V{x\ (5) 



is called k-th order supersymmetric quantum mechanics. In this approach there is a relation- 
ship between the supersymmetric 'Hamiltonian' H^s and the physical one H'^ = diag{/J, H} 
of polynomial type: 



H,, = l[{HP-e,). (6) 



1=1 



If one assumes that V{x) is a given solvable potential with normalized eigenfunctions i'n{x) 
and eigenvalues En, n = 0, 1, ... , equations (I4I6I) ensure that for any ipnix) such that 
B^il)n{x) 7^ it turns out that 

7 ( \ B^n{x) 

Wn\X) = —j^^^^^=^^^^^= (7) 

V [Bn — ei) . . . {En — Cfc) 

is a normalized eigenfunction of H with eigenvalue En- In general, the set {^„(x),n = 
0, 1, . . . } is not complete, since there can exist eigenstates i^eX^) of B. with eigenvalues 
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belonging as well to the kernel of B. By adding them to the previous set, the maximal set of 
eigenfunctions of H is thus given by: 

{V^,,(x),V^„(a;),i = = 0, 1,...} (8) 

The corresponding eigenvalues are {ej, En, i = I, k,n = 0,1, .. .}. 

In the maximal situation, the potential V{x) as well as the complete set of eigenfunctions 
of H are determined once the seed eigenfunctions Ui{x) of H (which not necessarily are 
physical) with eigenvalues e^, i = 1, . . . , k axe supplied. In particular, V{x) reads: 

V{x) = V{x) - {ln[W{u,, UkW, (9) 

W{ui, . . . , Uk) denoting the Wronskian of the seeds Ui{x), . . . , Uk{x). Let us illustrate the 
procedure more explicitly by means of the first and second order cases. 

2.1 First-order supersymmetric quantum mechanics 

Let us suppose that the intertwining operator is of first order 



d . . 
; h aiX) 



a/2 [ dx 

where the superpotential a{x) is to be determined. The use of equation @ leads to: 



(10) 



V{x) = V{x)-a'{x), (11) 
a\x) + a^{x)=2[V{x)-e], (12) 

i.e., a{x) must satisfy the Riccati equation (fT2l) . On the other hand, if a function u{x) such 
that a{x) = [\nu{x)]' is employed, equations (II 11121) become: 

V(x) = V{x) -[\nu{x)]", (13) 
Hu{x) = eu{x), (14) 

namely, u{x) obeys the initial stationary Schrodinger equation associated to e. 

Let us take now a solution a{x) (u{x)) to the Riccati (Schrodinger) equation (fT2l) ((fT4l)) 
for a fixed factorization energy e < Eq, where Eq is the ground state energy of H. Thus, 
equations (II 11131) indicate that the potential V{x) is determined completely, with a maximal 
set of normalized eigenfunctions {ijj^{x), ^/'„(x)} given by: 



ipt{x) oc exp 



a{y)dy 



ll)n{x) = (15) 



The corresponding eigenvalues are {e, En, n = 0,1, . . .}. Let us point out that the aim of 
the restriction e < Eq is to avoid that singularities appear in a{x), V{x) and also in the 
i/jeix), tpnix) of (fT5l) . Indeed, if e > £'o the seed solution u(x) will always have nodes in the 
x-domain of H and thus a{x) would have singularities at those points. If e < Eq, however, 
u{x) can have at most one zero. In particular, there is a subset of nodeless w-functions in the 
two-dimensional space of solutions associated to e < Eq, which will be used in the sequel 
for implementing the non-singular first-order SUSY transformations. 
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2.2 Second-order supersymmetric quantum mechanics 

Now, let the intertwining operator be of second order 

where ri{x), 7(x) are to be determined. Equation dH) leads to a set of equations relating 
V{x), V{x), r]{x), 7(0;) and their derivatives which, after some calculations reduce to: 

V = V-i, (17) 
I 2 

'y = ^ + ^-2V + d, (18) 

^ - Y + + T - 2^^' + c^r?' + c = 0, (19) 

with c, c? e M. For a given K(a;), the new potential V{x) and 7(x) are obtained from (|17ll8h 
once we find a solution 77 (x) of (fT9l) . which can be gotten from the Ansatz 

7]' = + 2f3r] + 2^. (20) 

By plugging (|20|) into (fT9l) , after some calculations we get = c and: 

P'{x) + P\x) = 2[V{x) - e], e={d + 0/2, (21) 

which is again a Riccati equation. We can work as well the related Schrodinger equation, 
which arises by substituting in (|2T)) = [lnM(a;)]': 

u" 

- — ^Vu = tu. (22) 

If c 7^ 0, ,^ takes the values i^c, and in this way we need to solve the Riccati equation 
(|2T)) for two factorization energies ei 2 = (rf ± ^/c)/2. Then one constructs algebraically 
a common solution 77(2;) of the corresponding pair of equations (|20|) . On the other hand, if 
c = one has to solve first the Riccati equation (|2TI) for e = d/2 and to find after the general 
solution of the Bernoulli equation resulting for ri{x) (see (|20|)). There is a clear difference 
between the situation with real factorization constants (c > 0) and the complex case (c < 0), 
suggesting to classify the solutions ri{x) based on the sign of c, which is next elaborated [30]. 

2.2.1 Real case (c > 0). 

Here we have ei 2 G M, ei 7^ 62, the corresponding Riccati solutions of (|2TI) being denoted 
by /3i,2(a;). The resulting formula for ri{x), expressed either in terms of /3i,2(a;) or of the 
corresponding Schrodinger seed solutions ^1,2(2^) becomes: 

^ 2(ei - 62) ^ 2(ei - £2)^1^2 ^ W^'(mi,»2) 

/3i(x)-/52(x) iy(wi,n2) Vr(wi,n2)' ^ ^ 

where (?) = fg' — is the Wronskian of / and g. It is clear from Eqs. (|17|23l ) that 
the new potential V{x) has no new singularities in (x^, x^) if W{ui,U2) is nodeless there. 
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The spectrum of H depends on weather or not its two 'mathematical' eigenfunctions -ipei 2 
associated to ei_2 which belong as well to the kernel of B can be normalized, namely 

Their explicit expressions in terms of Mi 2 are: 

V'ei oc — oc — -, ij,,oc — (x — -. (24) 

Ml W[Ui,U2) U2 W[Ui,U2) 

If both of them can be normalized, we arrive then to the maximal set of eigenfunctions of H: 

iJn = -J======= } . (25) 

Among the several spectral modifications which can be achieved through the real second- 
order SUSY QM, some cases are worth to be mentioned [15, 18]. 

(a) Deleting two neighbor levels. For 62 = -Ej-i, ci = Ei, U2 = ipi^i, ui = ipi, it turns 
out that the Wronskian is nodeless and V^^i, ^^2 non-normalizable. Thus, Sp(i/) = 
{£^0, ■ ■ • 5 -Ej-2, -Ej+i, . . . }, i.e., the levels Ei were 'deleted' for generating V{x). 

(b) Creating two new levels. For Ei_i < 62 < ei < E^, i = 1,2, . by taking U2, Ui with 
i + 1, i nodes respectively the Wronskian becomes nodeless, tp^^^, ipe2 normalizable and 
Sp{H) = {Eq, ■ ■ ■ , Ei_i, €2, ei, Ei,. . . }. 

(c) Isospectral transformations. They appear as a limit of the case in which two new levels 
are created for Ei^i < 62 < ci < Ei, when ui 2 satisfy either Mi,2(a^z) = or Ui^2{xr) = 0. 
In this case the Wronskian vanishes at xi or Xr, and ^^j, cease to be normalizable so that 
SpiH) = SpiH). 

2.2.2 Complex case (c < 0) [31]. 

Now e = ei G C, £2 = e, and since we look for V{x) real, it must be taken = Pi = 
P2{x). Hence, the real solution ri{x) of equation (fT9l) generated from the complex one (3{x) 
of (I2TI) becomes: 

, , 2Im(e) w'ix) , , W{u,u) 

Im[p(a;)J w[x) 2(e — e) 

Note that w{x) must be nodeless for x E {xi, Xr) to avoid new singularities in V{x). Since 
w{x) is non-decreasing monotonic {w'{x) = p), a sufficient condition ensuring the lack 
of zeros is 

lim u{x) = or lim u{x) = 0. (27) 
For transformation functions obeying (|27]) . V{x) is a real potential isospectral to V{x). 
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2.2.3 Confluent case (c = 0) [32,33]. 

We get now ^ = 0, e = ei = £2 G M; let us take a Riccati solution P{x) to (|2TI) for the given 
e. Thus, the general solution for the Bernoulli equation resulting of (|20l) reads: 



wix] 



Wo+ [ e'f^^'^^'^dx = Wo+ rHy)f dy, (29) 



where Xq is a fixed point in [x;, x^]. Once again, w{x) must be nodeless in order that V{x) 
has no singularities in Since w{x) is non-decreasing monotonic {w'{x) = [u{x)]'^), 

the simplest choice ensuring a nodeless w{x) is to take satisfying either 

pXQ 

lim m(x) = 0, /- = / [uiy)f dy < 00, (30) 

A, 

or 

pXr 

lim ■u(x) = 0, /+ = / [-ud/)]^ (iy < 00. (31) 

^^^-^ Jxo 

In both cases it is possible to find a Wo-domain for which w{x) is nodeless. The spectrum of 
H depends on the normalizability of the eigenfunction of H associated to e belonging as 
well to the kernel of B, with explicit expression given by: 

~ r]{x) u{x) 

iPe\X) OC — -— OC 



u{x) w{x) 

If it can be normalized, then the maximal set of eigenfunctions of H becomes: 

4),[x), %l)n[x) = — } . (32) 

Note that, for e > Eq, t ^ Em, m = 1, 2, . . . there exist solutions u satisfying (l30l) or (|3T1) 
such that 'ip^ is normalizable, i.e., the confluent second-order SUSY QM allows to embed a 
single level above the ground state of H. Moreover, since the physical eigenfunctions of H 
satisfy both (I30I31I) . they are also appropriate for implementing the confluent algorithm. Let 
us remark that, apparently, the first authors who realized that through the confluent SUSY 
QM it is possible to modify the excited part of the spectrum were Baye and collaborators 
[34, 35]. We thank one of the referees of this paper for this information. 

3 Trigonometric Poschl-Teller potentials and their SUSY 
partners 

Let us apply the previous techniques to the trigonometric Poschl-Teller potentials [16,27,29]: 

= + A..>1. (33) 

2sm (x) 2cos^(x) 



6 



Notice that, for 1/2 < A = z/ < 1, the V{x) of ([33]) is known as Scarf potential [9, 28]. The 
SUSY transformations for that periodic potential have been recently implemented [24] . 

Along the paper it will be extensively used the general solution of the Schrodinger equation 
Hu{x) = eu{x) for any positive value of the energy parameter e, which reads: 

u{x) = sin^(x) cos'' (x) I A 2^1 + ^/^, i - v^; A + sm^{x)] 

+Bsm^-^^{x) 2F1 + yf, i±|^ - v^; 1 - A; sm\x)] |, (34) 

where fi = X + u. We can find now the eigenfunctions ^„(x) of H, which satisfy the 
boundary conditions V^ri(O) = V^„(7r/2) = 0. Since ^/'n(O) = 0, it turns out that B = 
0. Moreover, for arbitrary e > the hypergeometric function involved in the remaining 
term diverges when x — 7r/2 stronger than the vanishing behavior induced by cos'^(x). In 
order to avoid this divergence so that ipn{7c/2) = 0, one of the two first parameters of the 
corresponding hypergeometric function has to be a negative integer, namely: 

u I E„ (u + 2nY 

|±^^ = -n ^-= 2 ' ^ = 0'1'2--- (35) 

By using the normalization condition it turns out that the eigenfunctions of H are: 



Mx) = J STraSi^i^ sin^(^) ^^^^i^) ^Fil-n, n+^^; A + i; sm'ix)]. (36) 



(i^+|)nr(A+|)r3(..+i) 

For implementing later the SUSY transformations, it is important to know the number 
of zeros of the Schrodinger seed solution which is going to be employed. These nodes 
depend on e, A, B (see expression (|34l) ). To determine that dependence, let us compare the 
asymptotic behavior of ^(a:;) for x 0, 7r/2. Indeed: 

m(x) ~ 5sin^-^(x), ~ (Aa + 56)cosi-'^(x), (37) 

„_ r(A+i)r(^^i) r(|-A)r(z.-i) 



r(t+yf)r(f-v^)' r(i+^+^)r(i+^-y|)- 

By asking that u{x) > when a; ~ 0, it turns out that B > 0. Without loosing generality 
let us take B = 1 and A = —b/a + q. Since for e < Eq u{x) just can have either one or 
zero nodes in (0, 7r/2), thus it will have one if g < while it will be nodeless if g > 0. For 
Eq < e < El, u{x) will have either two zeros for g < or just one for g > 0. In general, for 
Ei-i < e < Ei, u{x) will have either i + 1 nodes for g < or i ones for g > 0. 

Notice that the trigonometric Poschl-Teller potentials, and the corresponding Hamiltoni- 
ans, are invariant under the transformation x ^ n /2 — x, A — z/, z/ — > A. Its action onto 
the Schrodinger solution (|34l) . with a given e and specific values of the parameters (A, B), 
produces another solution with different parameters (Aai + BPi, Aa2 + BP2), where 

«i = - (M) b, «2 = a, 

^ r(i-A)r(|-,.) _ r(A-i)r(|-i.) 



r(i-t+v^)r(i-t-yf)' r(i+A^+^)r( 



This result will be used below to diminish the number of discussed SUSY transformations. 
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3.1 First-order SUSY partners 

Let us classify the first-order SUSY partners according to the changes induced on the initial 
spectrum. Three different cases have been identified [16]. 

(a) Deleting the initial ground state. Let us choose e = Eq and as seed the ground state 
eigenfunction of H, 

u{x) = ipo{x) oc sin'^(x) cos'^(x). (38) 
The SUSY partner potential of V{x) becomes: 

KW^^ + l^ll. A..>1. (39) 

2sm (x) 2cos"'(x) 

Since 'ipe{x) oc l/ipo{x) diverges at x = 0, 7r/2, the eigenvalues of H are given by (|35l) just 
with n = 1, 2, . . . , i.e., we have 'deleted' the ground state energy of H to generate V{x). 

The previous SUSY partner potential V{x) can be obtained of the initial one through the 
change v ^ u + 1, a property which is nowadays called shape invariance [9]. 

The fact that the singularities at a; = 0, 7r/2 are reinforced, increasing by one both parameters 
A, u, has to do with the vanishing at those points of the employed seed solution. This behavior 
is identical to the one observed at the origin for the singular term of the SUSY partners of 
effective radial potentials [6]. 

As an illustration, the potentials V{x) and V{x) for A = 3, = 4 are drawn in dashed 
and in gray respectively in figure 1 . 

(b) Creating a new ground state. Let us take now e < Eq and a nodeless seed solution u{x) 
given by (l34l) with B = 1, A = —h/a + g, g > 0. Since u{x) — > oo as x — > 0, 7r/2, then 
ipe{^) = V^e(7r/2) = 0, i.e., is a new eigenfunction of H with eigenvalue e. Note 
that Sp(if)={e, En, n = 0, 1, . . . }, namely, a new level has been 'created' at e for H. The 
singularities induced by on V(x) at x = 0, 7r/2 are managed by defining 

u{x) = sin"'^~'^(x) cos"'^~'^(x)f (x), (40) 

where v{x) is a nodeless bounded function in [0, 7r/2]. Thus we get: 

^(,),(^-2)(A-l)^(^-2)(.-l) ^^^^ 

^ ' 2sin2(x) 2cos2(x) i v ;j , , 

Notice that now the singularities at x = 0, 7r/2 are weakened, decreasing by one both param- 
eters A, z/. This is due to the divergence at both points of the employed seed solution, which 
once again is similar to the behavior at the origin for the singular term of the SUSY partners 
of effective radial potentials [3, 6, 10]. 

An example of the potential (f4T]) for A = 3, z/ = 4 is given by the black continuous curve 
of figure 1 . 

(c) Isospectral potentials. They appear from the transformations creating a new level at 
e < EqIW the limit when u{x) vanishes at one of the ends of the x-domain so that V'e(x) is 
not longer an eigenstate of H. In our example, two appropriate seeds are available, given 
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Figure 1: Trigonometric Poschl-Teller potential for A = 3, = 4 (gray curve) and its first order SUSY 
partners which arise from deleting the initial ground state Eo = 24.5 (dashed curve), creating a new ground 
state at e = 19 (black continuous curve), and making an isospectral transformation with the same e (dotted 
curve). 



by dUl) with A = l, B = OorA = -b/a, B = 1. In the first case m(0) = 0, and the 
corresponding divergence induced on V{x) can be handled by taking: 

u{x) = sin'^(x) cos^^'^(x)f (x), (42) 

v{x) being nodeless bounded in [0, vr/2]. With this choice it turns out that: 

Vix) = ^^+^-^^^^^-[lnvix)r, X>1, u>2. (43) 
2sm''(x) 2cos^(x) 

Since |^e(x)| oo when x ^ 0, then e ^ Sp(if ) and therefore H is isospectral to H. 

Notice the opposite changes of A, u suffered by the SUSY partner potentials V{x): the 
parameter A (u) is increased (decreased) by one since the seed solution vanishes (diverges) 
at a; = (a; = 7r/2). Once again this is similar to the modifications induced by SUSY on the 
term singular at the origin of effective radial potentials [3, 6, 10]. 

The potential (|43] ) for A = 3, z/ = 4 is illustrated by the dotted curve of figure 1. On 
the other hand, the second seed solution which satisfies u(n/2) = is obtained by changing 
x — ^ 7r/2 — x, \ u ^ A in (l42l). The corresponding isospectral SUSY partner potential 
arises from the same transformation applied to (|43l) . 



3.2 Second-order SUSY partners 

Let us explore the spectral modifications which can be induced in the three cases of the clas- 
sification of section 2 (a partial study is found in [29]). Our results suggest a rule which will 
be observed for the changes induced on the parameters A, u characterizing the singularities 
at X = 0, 7r/2 in the real and complex cases: if both seeds vanish (diverge) at x = 0, then 
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each one will increase (decrease) by one the parameter A so that at the end the coefficient of 
the divergent term of V{x) is obtained by making A ^ A + 2 (A ^ A — 2). On the other 
hand, if one solution vanishes while the other one diverges at a; = 0, then the correspond- 
ing singular term of V{x) will be the same as for V{x) (unchanged A). Something similar 
happens for the parameter u characterizing the singularity at a; = tt/2. This behavior is seen 
also for the singularity at the origin of the SUSY partners of effective radial potentials [6]. 



3.2.1 Real case. 

For ei 2 E M several possibilities of modifying Sp{H) are available. 

(a) Deleting two neighbor levels. Let us take ei = Ei, €2 = -Ej-i, ui{x) = V^i(x), U2{x) = 
ipi_i(x) (see equation (|36l)). It is straightforward to show that: 

W{ui, U2) oc sm^^+\x) cos^''+\x) W, (44) 

where 

^ _ W{2Fi[-i, z + /i; A + |; sin=^(x)], 2F1H + 1, i - 1 + /i; A + ^; sin^(x)]} ^^^^ 

sin(x) cos(x) 

is a nodeless bounded function in [0, vr/2]. The second-order SUSY partners of V{x) be- 
come: 

~ (A + l)(A + 2) (u + l)iu + 2) 
2sm (x) 2cos^(a;) 

The two mathematical eigenfunctions ip^^ oc U2/W{ui,U2), V'ea oc ui/W{ui,U2) of H 
associated to ei = Ei, £2 = -Ei-i do not obey anymore the boundary conditions to be 
physical eigenfunctions of H since 

Hm iV'ei ,(a;)| = 00. 

Thus, Sp(^) = {Eo, . . . E,_2, Ei+i, ...}. 

A plot of the potential (|46l) for A = 5, u = 8, generated by deleting the levels E2 = 
144.5, E3 = 180.5, is shown in dashed in figure 2, while the initial one is drawn in gray. 
Notice the stronger intensities of the singularities at a; = 0, 7r/2 of V{x) with respect to the 
corresponding ones of V{x) (compare the potentials (|33l) and (l46l)). 

(b) Creating two new levels. Let us choose now Ei_i < 62 < ei < Ei, and the corresponding 
seed solutions as given by (l34l) with B12 = 1, Ai^2 = — &i,2/'^i,2 + 11,2, 12 < 0, qi > 0, 
i.e., U2 and ui have i + I and i nodes respectively, making the Wronskian nodeless. In order 
to include the case when £2 < ei < Eq, let us assume that i = 0, 1, 2, ... , where we have 
introduced the formal fictitious level E_i = —00. It is important to 'isolate' the divergent 
behavior of the u solutions for x ^ and x — 7r/2 (see equation (l37l) ) by taking: 

ui^2{x) = sin^~^(a;) cos^~''(x)t;i,2(a;), (47) 
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vi,2{x) being bounded for x G [0,7r/2], f 1,2(0) 7^ 0, fi,2(7r/2) ^ 0. Since the second term 
in the Taylor series expansion of f 1,2(2^) is proportional to sin^(x), it turns out that v[ 2(0;) 
tend to zero as sin (a;) for x ^ and as cos (a;) for a; ^ 7r/2. A simple calculation leads to: 

W{ui, U2) = sin=^-2^(x) cos^-'^''{x) W, (48) 

where W = 1, i;2)/[sin(a;) cos(x)] is nodeless bounded in [0,7r/2]. The second-order 
SUSY partners of the Poschl-Teller potential (|33] ) are now: 

V-W.M^+(^Zi^_(,„wr A..>3. (49) 

2sm (x) 2cos''(x) 

Since 

lim ih^, „(x) = 0, 

then Sp{H) = {Eq, . . . , -Ej-i, €2, ei, Ei, . . .}, i.e., two new levels have been created between 
a pair of neighbor ones of H to generate V{x). 

A plot of the potentials (|49l ) for A = 5, u = 8, generated by creating the two new levels 
ei = 128, €2 = 115.52, is given by the black continuous curve of figure 2. Observe the 
weaker intensities of the singularities at x = 0, 7r/2 of V(x) compared with those of the 
initial potential (|33l) . 

(c) Isospectral transformations. They arise from those which create two new levels (see case 
(b)) in the limit when each seed vanishes at one of the ends of the x-domain. By simplicity, 
let us choose ui^2 as given in (|34l) with Bi^2 = 0, Ai^2 = 1 so that ^1,2(0) = 0. Since 
1^1,2(2;) I — > 00 when X 7r/2, it is convenient to express: 

''^1,2(3^) = sin'^(x) cos^^'^(x)fi,2(a;), (50) 

f i,2(x) being bounded in [0, vr/2], t'i,2(0) 7^ 0, f i_2(7r/2) 7^ 0. Once again, it turns out that: 

W{ui, U2) = sin^^+\x) cos=^-2'^(x) W, (51) 

where W = W{vi,V2)/[sm(x) cos(x)] is nodeless bounded in [0,7r/2]. The second-order 
SUSY partners of the Poschl-Teller potential are now: 

n.)^^2^^t^^-i^.^vr. A>i, .>3. ,52) 

2sm (x) 2cos^(x) 

Notice that 

lim I V'ei,2 (a;) I = 00, lim tfj,^ (x) = 0. 

This implies that ei 2 Sp(iJ), meaning that ^(x) is strictly isospectral to V{x). 

Note that a similar procedure for ui^2 satisfying Mi,2(7r/2) = can be applied. The cor- 
responding seed solutions and isospectral SUSY partner potentials are obtained by changing 
X 7r/2 — X, A — s> z/, z/ — > A in equations (|50II52I) . 

(d) Creating a new level. It appears from case (b) when one of the i + 1 nodes of U2 goes 
eitherto or to 7r/2. In the first case it is taken ,82 = 0, A2 = 1, Bi = 1, Ai = —bi/ai + qi. 
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qi > 0, so that ^2(0) = 0. In order to manage the singularity at a; = 7r/2 induced by ui^2 on 
V{x), it is convenient to write them as: 

Ui{x) = sin"'^~'^(x) cos"'^~'^(x)fi(a;), ^2(3;) = sin'^(x) cos^~^(x)f2(x), (53) 

'^'1,2(2^) being bounded in [0, tt/2], ^1,2(0) 7^ 0, fi^2(7r/2) 7^ 0. It can be shown that: 

W{ui, U2) = cos3-^''(x) W, (54) 

where W = W[sm^~^ {x)vi{x) , sin'^(a;)f2(x)]/ cos(a;) is nodeless bounded for x E [0, 7r/2]. 
The second-order SUSY partner potentials of V{x) are: 

I^W - 1^ On >Vr. A>1. .>3, (55) 

2sm (x) 2cos''(x) 

Since 

lim ipeiix) = hm 4'e2{x) = 0, hm |'?/'e2(a;)| = 00, 

thus Sp(iJ) = {Eq, . . . , Ei_i, ei,Ei, . . . }, i.e., we have embedded a new level ei in (Ei-i, Ei). 

The second possibility for generating a new level, in which U2{tt /2) = 0, can be obtained 
through the changes x ^ tt/2 — x, \ u, u ^ Xin formulae (I53H55I ). 

(e) Moving an arbitrary level. This can be achieved in the first place by taking the factoriza- 
tion energies as Ei_i = 62 < ei < E^ and the seeds in the way U2{x) = ^/',j_i(x), Ui{x) as 
given in (|34|) with Bi = 1, Ai = — &i/ai + gi, gi > so that ui{x) has i nodes in (0, vr/2). 
It is convenient to factorize the null and divergent behavior of the seed solutions Ui^2ix) at 
X = 0, 7r/2 by expressing them as: 

■ui(x) = sin^^'^(x) cos^^'^(x)f i(x), U2{x) = sin'^^x) cos'^{x)v2{x), (56) 

where 111,2(2;) are two bounded functions for x G [0, 7r/2], 111,2(0) 7^ 0, t'i,2(vr/2) 7^ 0. It 
turns out that W{ui, U2) is nodeless bounded for x G [0, 7r/2]. Moreover: 

lim ibe,{x) = 0, lim |'?/'eo(a;)| = 00, 

i.e., i^eiix) is an eigenfunction of H but V'e2(^) is not. The second-order SUSY partners of 
V{x) are given by: 

nx)=^:^^^,+^f^,-{HW{u,,U2W A,z.>l. (57) 
2sm (x) 2cos''(x) 

Since Sp(iJ) = {£"0, . . . , Ei_2, Ci, -Ej, • • • }, we conclude that the level Ei_i has been moved 
up to achieve ei. 

An example of the potentials (l57l) for A = 5, u = 8 is plotted in figure 2 (dotted curve). 
The initial level E2 = 144.5 has been moved up to achieve ei = 169.28. The 'intensities' of 
the singularities at x = 0, 7r/2 for V{x) remain the same as for the initial potential (|33l) . 
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Figure 2: Trigonometric Poschl-Teller potential for X — 5,iy — 8 (gray curve) and its second order SUSY 
partners (real case) which arise by deleting the levels E2 — 144.5, E3 = 180.5 (dashed curve), creating two 
new eigenvalues at ei = 128, £2 = 115.52 (black continuous curve), and moving the energy E2 — 144.5 up to 
ei = 169.28 (dotted curve). 



Another possibility is to take Ei^i < 62 < ei = Ei, the corresponding seed solutions in 
the way ui{x) = ipi{x), the U2{x) of (l34l) with A2 = — &2/«2 + 52, 92 < 0, i.e., ui{x) and 
U2{x) have i and i + 1 nodes respectively for x G (0, 7r/2). It is convenient to express 

= sin'^(x) cos'^(x)f i(x), U2{x) = sm^^'^(x) cos^^'^{x)v2{x), (58) 

^1^2(3;) being bounded for X G [0,7r/2], ^1,2(0) 7^ 0, fi,2(7r/2) 7^ 0. Once again, W{ui,U2) 
is nodeless bounded for x G [0, 7r/2]. Furthermore: 

lim lV'ei(a^)l = 00, Um il>e^(x) = 0, 

namely, ?/'e2 (x) is an eigenfunction of H while (2;) is not. The SUSY partner of V{x) is 
given as well by (l57l) . where now Sp(i/) = {Eq, . . . , -Ej-i, €2, -Ej+i, • • • }, meaning that the 
level Ei has been moved down to achieve €2- 

f/j Deleting an arbitrary level. This is attained of the previous case in the limit when the 
nonphysical seed acquires one zero at x = or s = 7r/2. For Ei_i = e2 < ei < Ei 
one possibility is to take U2{x) = Ui{x) as in (|34l) with Ai = 1, Bi = 0, so that 

ui{0) = 0. Thus 

'Ui(x) = sin'^(x) cos"'^~'^(x)f i(x), 'U2(a;) = sin'^(x) cos'^(a;)f2(a;), (59) 
f 1^2(2;) being bounded for x G [0, tt/2], t'i,2(0) 7^ 0, t'i^2(7r/2) 7^ 0. It turns out that: 

W{ui,U2) = sin2^+^(x) W, (60) 
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where W = [cos-'^ ''{x)vi{x), cos''(a;)f2(x)]/ sin(a;) is nodeless bounded for x G [0, vr/2]. 
Now we have 

lim ipei{x) = 0, Um |'?/'e^(a;)| = lim iV'eal^)! = ^5 

i.e., ei,2 ^ Sp(#) = {Eo, Ei_2, Ei, Ei+i, . ..}. The SUSY partner potentials of V{x) 
are given by 

?w.(A±iM,l!izi>:_(,„H.r. A..>i. (61) 

2sm (x) 2cos"'(x) 

It is seen that the level Ei^i has been deleted for generating V{x). 

Another option for deleting the level Ei^i can be achieved by changing x 7r/2 — x, A — > 
u, v ^ Xin equations 



3.2.2 Complex case. 

For e G C the solution u given in (|34|) is still valid, and the condition (|27|) required to avoid 
the zeros in the Wronskian can be accomplished in two ways. In the first place we make 
A = 1, B = and thus m(0) = while \u{x)\ ^ cx) as x — > n/2. The singularities induced 
on V{x) are handled by factorizing 

u{x) = sm'^{x) cos^^''{x)v{x). (62) 



Therefore: 



W{u, u) = sin'^+^(x) cos'-''^(x) W, (63) 

r.,.).(™±^,(^i^i)^_(,„W)". A>1, .>3. (64, 

2sm (x) 2 cos"' X 



W{v,v) 



2(e — e) sin(x) cos(x) 



The potentials V{x) of (|64l) and the Poschl-Teller initial one (|33l ) are isospectral. Their 
plots for A = 5, u = 8 are shown in figure 3, where the initial potential is drawn in gray 
while the dotted curve represents the one of (l64l) . 

Note that the second solution satisfying u{n/2) = 0, lim^-^o 1^(2^)1 ~^ 00, and the cor- 
responding SUSY partner potential V{x), are obtained by changing x n/2 — x, A — > 
u, 1/ — i> A in equations 



3.2.3 Confluent case. 

For e = ei = £2, several possibilities of modifying the initial spectrum appear. 

(a) Creating a new level. Let us choose M 3 e ^ Ei, for which two seeds are available for 
implementing the confluent algorithm. The first one arises by taking A = 1, _B = in (|34l) : 

u(x) =sin^(x) cos''(x) 2F1 (f + A/f, a/I; A + i; sin2(x)) =sin^(x) cos^"''(x)t;(x) (65) 
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Figure 3: Trigonometric Poschl-Teller potential for A = 5, = 8 (gray curve) and its second order SUSY 
partner (complex case) which arises by using e = 176.344 + 1.5i with a seed vanishing at the origin (dotted 
curve). 



v{x) being bounded for x E [0, tt/2], v{0) ^ 0, v{tx /2) ^ 0. The calculation of the integral 
of equation (|29l) with xq = leads to: 



W[X) - Wo + 2. (A+i)„m!(2A+2m+l) 



X 



.i^2(i 



m=0 



A+m+|; A + |, A+m+|; sin^(x)) . 



2' 



2' ' 



(66) 



Notice that w{x) is nodeless in [0, 7r/2] for > while it will have one node for wq < 0. 
Let us choose a nodeless w{x), as given in (l66l) with u'o > 0. Its divergent behavior for 
X 7r/2, being of kind cos^~^'^(x), will change the coefficient of the second term of the 
Poschl-Teller potential (1331 ). so it is convenient to factorize 



w{x) 



cos 



3-2!/, 



x) W(x) 



(67) 



W(x) being nodeless bounded for x E [0, 7r/2]. The confluent second-order SUSY partner 
potentials of V{x) become: 



V{x) 



(A-l)A (z/-3)(z/-2) 



2 sin ix) 



2 cos^fx) 



{ln[W(x)]}", A>1, z/>3. 



Since il)e{x) oc u{x)/w{x) satisfies: 



lim ibeix) 



0, 



(68) 



(69) 



then Sp(i/) = {e, n = 0, 1, . . . }, e ^ 

As an illustration, in figure 4 we have drawn in gray the initial potential for A = 5, z/ = 8 
and its SUSY partner (|68] ) by the black continuous curve. It is seen the different intensities 
of the singularities for both potentials at x = 7r/2. 
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Figure 4: Trigonometric Poschl-Teller potential for A = 5, = 8 (gray curve) and its second order SUSY 
partners (confluent case) which arise from creating a new level at e = 147.92 (black continuous curve), making 
an isospectral transformation with e = 162 (dotted curve) and deleting the eigenvalue — 180.5 (dashed 
curve). 



Notice that the second seed, which is appropriate to implement the confluent algorithm, 
and the corresponding SUSY partner potential, are obtained by changing x ^ 71/2 — x, A— > 
u, ^ A in equations (I65II68I) . 

(b) Isospectral transformations. They appear in several different ways, in the first place as 
two limits of the previous case when the eigenfunction of H associated to e ceases to satisfy 
the right boundary conditions. This happens, e.g., if we take u{x) as in (l65l) and the w(x) 
of (|661 ) with wq = 0. Besides the divergent behavior of w{x) as x ^ n/2, it turns out that 

w{x) ^ as sin^'^+^(a;) when a; — > 0, so that: 

w{x) = sin2^+^(x)cos3-2''(x) >V(x), (70) 
being nodeless bounded for x E [0,n/2]. The SUSY partner potential of V{x) is: 

2sm (x) zcos'^^x) 



Notice that: 



lim = oo, lim ipe{x) = 0, (72) 

a;— >0 x— >5 



i.e., e ^ Sp(if) and therefore H has the same spectrum as H. 

An example of the potentials (TtTI) for A = 5, u = 8, e = 162 is shown in dotted in figure 
4. It can be seen that the stronger intensity of the singularity at x = of ^(2;), compared 
with V{x), is 'compensated' by its lower value at x = n/2. 

A second alternative to produce isospectral potentials consists in changing x 7r/2 — 
X, A — > z/, z/ ^ A in equations (I65H68I) and taking i^o = in the resulting formulas. The 
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corresponding SUSY partner potential is obtained by substituting a; 7i/2—x, X u, z/ — 
A in equations (|70I71I) . 

The third confluent isospectral transformation uses as seed a physical eigenfunctions of 
H, i.e., e = Ei, u{x) = -ipiix). The expression for w(x) is obtained from (|66l ) by realizing 
that the solution (|65] ) is proportional to the eigenfunction (|36l ) when e ^ E^, 

ipi = Q Jhn sin^(x)cos''(x) 2E1 (f + y |' ^ " V ^ ^ h ^^^^ ^^^J ' 



C,; 



2{fi + 2i)i\T{fx + i){\ + ^)i 



Moreover, in this limit the infinite summation of (l66l) truncates at m = z, so that: 

nj{x) =Wo + C\ E ^^A+iv!n!(2A+2.n+l) 
m=0 ^ ^' 

X 3i^2 (I - - I + ^ + ^' ^ + + I' ^ + I' ^ + + I' sin^(a;)) • (74) 

If Wo > or Wo < — 1, w(x) is nodeless bounded for x G [0, 7r/2]. Now there is not change 
in the intensities of the singularities at x = 0, 7r/2 for 'V[x\ namely: 

y^^) = 7^^^, + - {ln[^(x)]}", A, > 1. (75) 

2sm ix) icos^ X 



It turns out that: 



lim ^pJx) = 0, (76) 



i.e., e = Ei e Sp(if) and thus H and H axe isospectral. 

(c) Deleting an arbitrary level. This case appears in the limits as wo 0, —1 of the 
isospectral transformations involving as seed the physical eigenfunction ^t(x). For wq — > 0, 
w{x) ~ sin^'^^"^(x) when x ^ so that: 

w{x) = sin2^+^(x) W(x), (77) 

where >V(x) is nodeless bounded in [0, 7r/2]. Since 

lim |'?/'e(x)| = cx), \imijj^{x) = 0, (78) 

then e = Ei^ Sp(^) = {Eq, . . . , Ei+i, ...}. The SUSY partner potential of V{x) is: 

Vix) = (^ + ^^1) + '^ + - {ln[W(x)]}", A, . > 1. (79) 

^ ^ 2sin2(x) 2cos2(x) ^ l v ;j; > > 

It is seen that we have deleted the level Ei to produce V{x). 

An illustration of the potentials (|79l) for A = 5, z/ = 8 is shown in dashed in figure 4. The 
deleted level is E-^ = 180.5, and the intensities of V{x) and V (x) at x = differ as predicted 
by equations ([33]) and (l79l) . 

The case when wq —1, which leads also to the deletion of the level E^, can be achieved 
from equations (I77I79I ) by the change x ^ n/2 — x, \ ^ u ^ X. 
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4 Conclusions 



The supersymmetric quantum mechanics of first and second order have been used to gener- 
ate new exactly solvable Hamiltonians departing from the trigonometric Poschl-Teller poten- 
tials. Several interesting possibilities to modify the initial spectrum have been studied, and 
it has been shown that the deformations induced by the second order algorithm can be non 
standard, in the sense that the main spectral changes appear above the ground state energy of 
the initial Hamiltonian. Specifically, we have shown that a pair of levels can be embedded 
between two neighbor initial ones. It has been possible also to delete two neighbor energies. 
Specially interesting is the possibility of embedding a single level at any arbitrary position. In 
addition, it is possible to move up or down a generic physical energy as well as to delete it. It 
is worth to notice that some spectral modification can be achieved in several different ways. 
For example, the strictly isospectral mappings can be obtained through the real, complex 
and confluent second-order transformations (see the potentials in (I52I64I7 11751 )). However, 
if we want to produce an isospectral potential such that the coefficients of the singularities at 
X = 0, 7r/2 would be changed in a specific way, then the number of options becomes smaller. 
In particular, if the isospectral SUSY transformation is not going to modify the intensities of 
the two singularities at a; = 0, 7r/2, then we will have to apply a confluent transformation in- 
volving as seed a physical eigenfunction of the trigonometric Poschl-Teller Hamiltonian (see 
equation (l75l)). A similar discussion could be elaborated for the other cases having several 
possibilities to achieve the same final spectrum. Our general conclusion is that the super- 
symmetric quantum mechanics is a powerful mathematical tool for designing potentials with 
an arbitrarily prescribed spectrum. 
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